Quantum versus classical phase-locking transition in a driven- chirped oscillator 
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Classical and quantum- mechanical phase locking transition in a nonlinear oscillator driven by a 
chirped frequency perturbation is discussed. Different limits are analyzed in terms of the dimen- 
sionless parameters Pi = e/ \/2mhjjQa and P2 = (3/i/3)/(4mv^) (e, a, /3 and ujq being the driving 
amplitude, the frequency chirp rate, the nonlinearity parameter and the linear frequency of the 
oscillator). It is shown that for P2 <^ Pi + 1, the passage through the linear resonance for Pi above 
a threshold yields classical autoresonance (AR) in the system, even when starting in a quantum 
ground state. In contrast, for P2 ^ Pi + 1, the transition involves quantum-mechanical energy 
ladder climbing (LC). The threshold for the phase- locking transition and its width in Pi in both AR 
and LC limits are calculated. The theoretical results are tested by solving the Schrodinger equation 
in the energy basis and illustrated via the Wigner function in phase space. 
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I. INTRODUCTION 

Autoresonance (AR) is a generic nonlinear phase- 
locking phenomenon in classical dynamics. It yields a 
robust approach to excitation and control of nonlinear os- 
cillatory systems by a continuous self- adjustment of sys- 
tems' parameters to maintain the resonance with chirped 
frequency perturbations. Applications of AR exist in 
many fields of physics, examples beingatomic and molec- 
ular systems [ll, nonlinear optics [sl, Josephson junc- 
tions [4], hydrodynamics fs*], plasmas [6^, nonlinear waves 
[7], and quantum wells [8]. Most recently, AR served as 
an essential element in the formation of trapped anti- 
hydrogen atoms at CERN [l^ and in studying the 
effect of fiuctuations in driven Josephson junctions prj . 
While the classical AR is well understood, the investiga- 
tion of the quantum-mechanical limits of the problem has 
started only recently [l, |Tl|, |l2| • The present study fo- 
cuses on the interrelation between the classical and quan- 
tum descriptions of the autoresonant transition in the 
simplest case of a driven Duffing oscillator (modeling a 
driven diatomic molecule [13] or a Josephson junction 
for example) governed by the Hamiltonian 
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where cpd = J ojddt^ ujd = oJo — od is the chirped driving 
frequency and a, /3 > 0. We wih assume that initially 
our oscillator is in a thermal equilibrium with the en- 
vironment at temperature T, but the chirped system's 
response is sufficiently fast to neglect the effect of the 
environment on the out-of-equilibrium dynamics pj|. 

Classically, in autoresonance, after passage through 
the linear resonance at t = 0, the driven oscillator grad- 
ually self- adjusts its oscillation frequency to that of the 
drive by continuously increasing its energy yielding 
a convenient control of the dynamics by variation of an 



external parameter (the driving frequency). The transi- 
tion to the classical AR by passage through linear reso- 
nance has a threshold on the driving amplitude, scaling 
as ~ [14]. This threshold is sharp if the os- 

cihator starts in its zero equilibrium, but in the presence 
of thermal noise it develops a width, scaling as T^/^ [l^. 
Both the AR threshold and its width have their quantum- 
mechanical counterparts, which will be discussed in this 
work. 

When the problem of autoresonant transition is dealt 
with quantum-mechanically, two questions must be ad- 
dressed. First, what are the differences between the clas- 
sical and quantum evolutions of the chirped-driven non- 
linear oscillator? In dealing with this question, Ref. [l2| 
suggested that the natural quantum-mechanical limit of 
the classical AR is a series of successive Landau-Zener 
(LZ) transitions or energy ladder climbing (LC), 
where only two adjacent energy levels of the driven oscil- 
lator are coupled at any given time. In contrast, the clas- 
sical AR behavior takes place when many levels are cou- 
pled at all times during the excitation [l7j . We will adopt 
and further develop this point of view here and describe 
different regimes in the problem in terms of two dimen- 
sionless parameters Pi, 2 suggested in [12]. These param- 
eters are defined via the three physical time-scales in the 
system, i.e. the inverse Rabi frequency Tr = \/2mhjjjQ/e^ 
the frequency sweep time scale = 1/ ^/a^ and the char- 
acteristic nonlinearity time scale T^l = {3h/3) / {Ama) 
(the time of passage through the nonlinear frequency 
shift between the first two transitions on the energy lad- 
der). Then, by definition 
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(a measure of the nonlinearity in the problem). We will 
show in this work that this parameter space describes all 
limiting cases of quantum-mechanical evolution in our 
system, including quantum initial conditions, the subse- 
quent transition to either LC or AR, and the associated 
threshold phenomenon. Note, that Pi, 2 have a meaning 
only in the case of a chirped system, because of the new 
time scale, T^, associated with this case. 

The second question, which must be addressed in the 
quantum-mechanical formulation of our problem is that 
of quantum fluctuations. As mentioned above, in the 
presence of thermal noise, the classical AR transition 
probability develops a width, scaling as T^/^ with tem- 
perature [15]. Nevertheless, at very low temperatures, 
the quantum fluctuations should be taken into account. 
Recent experiments by Kater et al. [11] demonstrated 
quantum saturation of the width of the phase-locking 
transition in superconducting Josephson junctions at suf- 
ficiently low temperatures, confirming the prediction that 
T in the classical width formula should be replaced 
by an effective temperature, Teff, where Teff = T for high 
temperatures and saturates at Teff = huJo/2kB at low 
temperatures. The experimental results imply that the 
fluctuations only determine the initial conditions of such 
a non-equilibrium oscillator and do not affect its time 
evolution. In this work, we will address the effect of 
quantum fluctuations in the AR problem theoretically 
and provide further justification of using the classical AR 
threshold width formula with T replaced by Teff. 

The scope of the paper will be as follows. In Sec. 
II we will use the quantum-mechanical energy basis in 
the rotating wave approximation and compare the driven 
dynamics of our oscillator in the quantum and classical 
regimes numerically. Section III will present the analytic 
description of the transition to phase-locking in terms of 
the Pi, 2 parameter space in both classical AR and quan- 
tum LC regimes. In the same Section, the theory will 
be compared with numerical simulations. Section IV will 
focus on the effect of quantum fluctuations on the width 
of the phase-locking transition. Finally, we will address 
the phase space dynamics in the problem in Sec. V by 
solving the quantum Liouville equation for the Wigner 
function numerically and compare the phase space evo- 
lution with that in the energy basis. Our conclusions will 
be summarized in section VI. 



II. CHIRPED DYNAMICS IN THE ENERGY 
BASIS 

We write the wave function of the oscillator governed 
by Eq.Q, = Xln^nl^n), in the energy basis |V^n) of 
the undriven {e = 0) Hamiltonian ([1]). The associated 
Schrodinger equation yields 



where we approximate the energy levels [18| 

En ^ hujo (n + 1/2 - /^^(n^ + n + 1/2)) , (5) 

n = 0, 1, 2, (3q = and e = e^f^^. We assume a 

weak coupling, £ <C £^0, and, consequently, neglect the 
nonhnear correction of order e/3q/huJo in the coupling 
term. Next, we define Cn = e^'^'^^Cn, where ujn = En/h, 
substitute this definition into Eq. (|4]), and neglect the 
nonresonant terms (rotating wave approximation) to get 

dt 2^2^"^ 

+V^Cn-ie^^--^-^-^^^), (6) 

where uJn,n+i = ^n+i - uJn = ^0 - 2ujo{n + l)/3q. Fi- 
nally, we introduce Bn = (7^6"*-^^'^^^, where 7^ = 
nat — n{n -\- l)uJo/3q and the dimensionless slow time 
r = ^/at^ associated with the change r^/2 of the driving 
phase due to the driving frequency chirp. Then Eq. (j6j) 
can be written as 

= Tn^n + Y + V^Pn-l) , (7) 

where = n[r — (n + l)P2/2], and Pi = e/^s/2ahuJom^ 
P2 = 2bjQ^^ql ^fa^ as defined in the Introduction. Note, 
that Pi characterizes the strength of the coupling be- 
tween the adjacent levels, while P2 is associated with the 
nonlinearity in the problem and determines the degree of 
classicality in the system (see Sec. V). Note also that the 
rotating frame here is chirked instead of the usual, fixed 
frequency frame and, thus, there remains an explicit time 
dependence in Eq. ([7]). Our goal is to analyze these slow 
evolution equations, but first, we discuss different limits 
in the driven system in Pi ^2 parameter space. 

The comparison between the classical AR and the 
quantum LC regimes was first discussed by Marcus et al. 
[12], who suggested the nonlinear resonance classicality 
criterion, P2 <C Pi, by requiring that the classical reso- 
nance width would include more than two quantum lev- 
els. Since the chirp rate cancels from this creterion, the 
latter characterizes the nonlinear resonance phenomenon 
in the system driven by constant frequency drive as well. 
The chirping introduces a new effect, i.e. a possibility of 
a continuous self- adjustment of the energy of the oscilla- 
tor to stay in resonance with the drive. This yields a new 
condition, separating the classical AR and quantum LC 
transitions, where the dynamics of the chirped system is 
very different. In the LC transition, only two levels are 
coupled at a time and the system's wave function climbs 
the energy ladder by successive LZ transitions [iG^. For 
example, Eq. ([7]) yields the following two- level transfor- 
mation matrix for the n — 1 ^ n transition 

/ (n - l)r - -^P^ ^P, \ 
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FIG. 1: (color online) The dynamics in the energy basis (a-d) and the corresponding phase space dynamics of the Wigner 
function (e-h) in the quantum ladder climbing regime, P2 = 8. The subplots correspond to times r = (a,e), 30 (b,f), 60 (c,g), 
and 90 (d,h). Only a single level is highly populated in the phase- locked group of levels. The dashed lines in subplots (e-h) are 
the separatrices of the external potential well. The dimensionless phase space coordinates are rescaled as ^ = \/px, v — \fp)U. 



We can calculate the time of the n^^ transition, 
by equating the diagonal elements in this matrix, i.e. 
Tn = nP2^ SO the time interval between two successive 
transitions is Ar = P2- On the other hand, the typi- 
cal duration Ar^z of each LZ transition has two distinct 
limits [19,]. In the non adiabatic (sudden) limit (Pi <C 1), 
Atlz is of the order of unity, while in the opposite (adi- 
abatic) limit, Atlz ^ Pi- Therefore, by comparing Ar 
and Atlz^ we expect to see well separated successive LZ 
steps, i.e. the LC, provided P2 ^ ^1 + 1 which describes 
both the sudden and the adiabatic limits. In contrast, the 
classical AR transition requires P2 <C Pi + 1, which co- 
incides with the nonlinear resonance classicality criterion 
mentioned above, when Pi 1. In section V, a differ- 
ent argument will be suggested to explain why classical 
mechanics yields the correct description of the transition 
to autoresonance when a stronger inequality, P2 <C 1, is 
satisfied, even when the system starts in the quantum 
mechanical ground state. Next, we discuss the numerical 
solutions of the problem and compare different regimes 
of chirped-driven dynamics. 

We have solved Eqs. ([7]) numerically, subject to ground 
state initial conditions Bn{ro) = (5n,o at tq = — 8 (the lin- 
ear resonance corresponds to r = 0). Each of the Figs. 
[IH3] corresponds to a different value of the nonlinearity 
parameter P2 and show the distribution of the popula- 
tion of the levels in the system at four different times 
(subplots a-d). The subplots e-h in the Figures show the 
associated Wigner distributions (see Sec. V) at the same 



times. Figure [T] shows the case of the LC dynamics for 
P2 = 8 and Pi = 0.8 at r = 0, 30, 60, and 90 (subplots 
a-d), and illustrates a clear time separation beyond the 
linear resonance between the successive LZ transitions. 
For example, we observe two groups of resonant and non- 
resonant levels at r = 90, separated by a valley centered 
at about n = 6. We find that the location of the resonant 
levels is determined by the slow time, i.e. n ^ r/P2, as 
shown above. Thus, the resonant (phase- locked) state 
in the system is efficiently controlled via the driving fre- 
quency and a given final state can be reached (and main- 
tained) by terminating the frequency chirp at the desired 
energy level. We also see that there exists a single highly 
occupied level in the resonant group of levels at any given 
time, indicating successive LZ transitions, as expected in 
the LC regime. 

Our second numerical example is presented in Fig. 
[2] and illustrates the intermediate regime (as discussed 
above) with Pi = P2 = 1 and r = 0, 8, 16, and 24 (the 
subplots a-d). As in Fig. (TJ a clear separation between 
the resonant (n < 5) and nonresonant (n > 20) groups 
of levels is seen in the Figure. We see that, typically, 
several levels are excited in the resonant group, but their 
number is small, so the driven dynamics can not be con- 
sidered as classical. The last example (see Fig. [3]) cor- 
responds to the classical regime, P2 = 0.2, Pi = 1.9 and 
r = 0,4,8,12. One observes a separation between res- 
onant and nonresonant groups at r = 12. Note that in 
all our numerical examples about 50% only of the initial 
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FIG. 2: (color online) The dynamics in the energy basis (a-d) and the corresponding phase space dynamics of the Wigner 
function (e-h) in the intermediate regime, P2 = 1. The subplots correspond to times r = (a,e), 8 (b,f), 16 (c,g), and 24 (d,h). 
Few levels are simultaneously excited in the phase-locked group in the intermediate regime. The dashed lines in (e-h) are the 
separatrices of the external potential well. The dimensionless phase space coordinates are rescaled as ^ = \/~^x^ v — \/pu. 
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FIG. 3: (color online) The dynamics in the energy basis (a-d) and the corresponding phase space dynamics of the Wigner 
function (e-h) in the classical AR regime, P2 = 0.2. The subplots correspond to times r = (a,e), 4 (b,f), 8 (c,g), and 12 (d,h). 
Many levels are excited in the phase- locked group. The dashed lines in (e-h) are the separatrices of the external potential well. 
The dimensionless phase space coordinates are rescaled as ^ = \/px, v — \fpu. 
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state is transferred to the continuing phase-locked state, 
leading to the question of resonant capture probability, 
which is discussed next. 



III. RESONANT CAPTURE PROBABILITY 
A. Threshold for phase-locking transitions 

For a given set (Pi, P2) we define the resonant capture 
probability, 

00 

P=Y.\Bn? (9) 

n=nc 

where Uc is the number of the level separating the reso- 
nant and nonresonant groups of levels at sufficiently large 
times. For a given value of P2, the probability P depends 
on the driving parameter. Pi. For example, in the case in 
Fig. [U we use nc = 6 and the resonant capture probabil- 
ity is P = 0.48. Similarly, in the two examples in Figs. 
2 and 3, we choose = 10, 40 to get P = 0.62, 0.66, 
respectively. 

We calculate the resonant capture probability by solv- 
ing Eqs. ([7j) numerically subject to initial conditions, 
Bni^o) = ^n,o (the ground state), for different values of 
Pi ^2 and To = —10. For a fixed P2, the capture probabil- 
ity P is a monotonically increasing, smoothed step func- 
tion of Pi. We define the threshold for efficient phase- 
locking transition, Pf^, as the value of Pi for 1/2 cap- 
ture probability, i.e. P{P^^) = 0.5. The full circles in 
Fig. Hlshow Pf^ for different values of P2. The dashed 
and dashed-dotted lines are the assymptotic theoretical 
predictions for the quantum LC and classical AR (see be- 
low), which agree with the results of our simulations in 
both limits. The line P2 = Pi + 1 is the separator between 
the classical and the quantum regimes of the chirped non- 
linear resonance, as discussed in Sec. II. This line crosses 
the threshold line Pf^ at (Pi,P2) ~ (0.8,1.8). One can 
see in the Figure that indeed, this point separates very 
different dependences of Pf^ on P2 associated with the 
quantum and classical dynamics of the chirped system. 
One can also see the oscillating pattern of the threshold 
Pf^ at 1 < P2 < 5, where the transition to phase- locking 
involves a mixture of LC and multi-level LZ steps. Next, 
we calculate the threshold for phase-locking transitions 
analytically. 
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FIG. 4: (color online) Different regimes of phase-locking tran- 
sition in the chirped oscillator. The full circles show the lo- 
cation of the numerical 1/2 resonant capture probability (the 
threshold for phase-locking transition) obtained by solving 
the Schrodinger Eq. ([Tj) subject to initial condition in the 
ground state. The dashed and dashed-dotted lines represent 
the theoretical thresholds in the LC and AR regimes, respec- 
tively. The line P2 = Pi + 1 separates the classical AR and 
the quantum LC regimes. 

this case, we can calculate the probability of each transi- 
tion separately, and multiply the probabilities. The two 
level transformation matrix (|8]) in the energy basis for 
the n — l^n transition yields the transition probability 
via the LZ formula 

Pn-l^n = 1 - r", (10) 

where r = e~^^i . We define the probability P for cap- 
ture into resonance in this case as the probability of occu- 
pying a sufficiently high energy level N after N successive 
LZ transitions, i.e. 

N 

P=l[{l-r% (11) 

Then, solving P(r) = 0.5, one finds the threshold for the 
LC transition. 



B. Quantum- mechanical ladder climbing 

In the quantum LC regime the nonlinearity parameter 
P2 determines the time interval between successive reso- 
nances [see Eq. (|8])]. In the case of a strong nonlinearity, 
at any given time only two levels are coupled, and the dy- 
namics can be modeled by successive LZ transitions. In 



Pf^ = 0.79, (12) 

where for two digits accuracy we used = 5 in the 
rapidly converging product (pT]) . Thus, the capture into 
resonance occurs in the first few LZ transitions and one 
can choose nc = 5 (see Fig. [T]in the definition Eq. ([9]) 
for calculating the capture probability near the threshold. 
This prediction is valid for large P2, as mentioned above. 
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The dashed Hne in Fig. |4] represents Eq. (p!2]) . while the 
numerical result for 1/2 capture probability is shown by 
full circles. One can see a very good agreement between 
the two results in the LC limit (P2 > 5). However, in 
the intermediate regime (1 < P2 < 5), oscillations in Pf^ 
are observed before convergence at the predicted LC Hne. 
These oscillations are due to the mixing of more than two 
neighboring levels in passage through resonance (see Fig. 

Ej). 



C. Classical autoresonance 
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As P2 decreases, a growing number of levels are cou- 
pled simultaneously and the dynamics becomes increas- 
ingly classical. The classical AR phenomenon is now 
well understood [14]. If one starts in the zero amplitude 
equilibrium, the autoresonant phase-locking is achieved 
for drives of amplitude e above the critical value Scr = 
1.34a^/^/3~^/^mco'y^ [14]. When expressed in terms of 
Pi 2, this classical threshold is translated into 



Pf^ = 0.82/V^. 



(13) 



When thermal fluctuations are included, the transition 
probability develops a width scaling as T^/^ with tem- 
perature [15|. At the same time, the threshold for 1/2 
capture probability remains the same. Thus, Pf^ in 
Eq. (p!3|) is the classical counterpart of the quantum- 
mechanical observable P{^^ in Eq. (fT2]) . This classical 
threshold is shown in Fig. |4]by dashed-dotted line, illus- 
trating excellent agreement with simulations (full circles) 
in the classical regime, P2 <C 1. It should be emphasized 
that the simulation results in the Figure are solutions of 
the quantum- mechanical equations ^ with parameters 
in the classical regime, while the probabilities of capture 
were calculated using the proper transition level Uc for 
each value of P2, as defined in Eq. (j9j). In the next Sec- 
tion, we discuss the width of the autoresonant transition. 



IV. THE WIDTH OF THE PHASE-LOCKING 
TRANSITION 

Another observable of the phase-locking transition 
mentioned above is the width of the transition, which 
we define as the inverse slope {dP/dPi)~^ of the phase- 
locking probability at P = 1/2. This width depends 
on the initial conditions governed by the thermal equi- 
librium with the environment. Classically, the thermal 
width of the autoresonant transition scales as ^] 



As = l.2?>^amkBT. 



(14) 



However, at very low temperatures, the classical ther- 
mal noise becomes negligible, but quantum fiuct nations 
remain. Recent experiments in Josephson circuits prj 



FIG. 5: (color online) The width of the phase-locking tran- 
sition in passage through resonance. Numerical results (full 
circles) are compared with the theoretical predictions (solid 
lines) in the classical AR (on the left) and the quantum LC 
(on the right) regimes. The system was in the ground state 
initially. 



demonstrated quantum saturation of the transition width 
at the value obtained from Eq. (p!4|) . but with ksT re- 
placed by the energy of the ground level. More gen- 
erally, it was suggested to calculate the width by replac- 
ing T in the classical formula by an effective temperature, 
-^eff ~ It^ CQth( ^^^°y ), in agreement with the experimen- 



tal results. Using Teff, we can translate Eq. ([11]) into the 
transition width in terms of Pi 



APf = 1.23V^BTeff/2ncJo, 



yielding 



APf 



0.61 



(15) 



(16) 



in the zero temperature limit. The Josephson circuit ex- 
periments [11] were performed with P2 = 0.00053, i.e. 
well inside the classical region (see Fig. |4]). Interest- 
ingly, these experiments allowed to characterize the ini- 
tial quantum "temperature" Teff of the system by mea- 
suring the final classical autoresonant state of the chirped 
excitation. We will justify this approach in the next Sec- 
tion by analyzing the dynamics of the associated Wigner 
function in phase space. In contrast to Eq. (p!5|) valid 
when the final state of the system is classical (P2 <C 1), 
the threshold width of the phase-locking transition in the 
LC regime (P2 ^ 1) can be calculated by evaluating the 
slope of P(Pi) from Eq. ([Ilj) at Pi = Pf^ = 0.79, yield- 
ing 



APf 



0.66, 



(17) 



where we assume that the system is in the ground state 
initially. Figure [5] summarizes our theoretical predictions 
for the width of the phase- locking transition (for the same 
parameters as in Fig. [4j) and compares them with those 
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from numerical simulations via the Schrodinger equation 
([7j). We see a good agreement in both the AR and LC lim- 
its, but notice significant oscillations of the width in the 
intermediate range of P2. Remarkably, while the thresh- 
olds in the classical and quantum-mechanical limits have 
very different scalings, the widths of the transitions are 
nearly the same. 



V. CHIRPED DYNAMICS IN PHASE SPACE 

Phase space dynamics comprises a convenient frame- 
work for comparison between classical and quantum evo- 
lution of the system. The Wigner function is one of the 
most useful phase space representations of the quantum 
mechanics, since it reduces to the classical phase space 
distribution in the limit of ^ 0. In this Section, we 
will study the dynamics of the Wigner function in our 
chirped oscillator problem in both the fixed and the ro- 
tating frames and discuss the transition to the classical 
limit in the problem. 



A. Wigner dynamics in the fixed frame 

The Wigner function associated with the ID 

2 

Hamiltonian of form H{x^p) = ^ + V{x^t) is governed 
by the quantum Liouville equation 

df , S/ 1 dVdf d^'+^Vd^'+^f 



zAj^u—-—-—— — 



dt ' ^ dx m dx du ^ m(2/ + l)! dx^^^^ dv?^^^'' 

(18) 

where u = p/m and we neglect possible decay and de- 
coherence processes. We take a low temperature limit, 
neglect the nonlinearity initially, and assume that the 
initial state of the system is in equilibrium with the en- 
vironment, i.e. [2Q|. 



mujo 

27r/CBPeff 



(19) 



where Teff = (fia;o/2/c^) coth(ficjo/2/cB^) is the effective 
temperature. Note that Teff ^ T at high temperatures, 
while Teff huJo/2kB at T ^ 0. 

In the case of interest the potential is a quartic [see 
Eq. ([1])] and, therefore, only one term survives in the 
right hand side of (p!8)) . allowing to rewrite this equation 
in the following dimensionless form 



and € = e/mLujQ. In addition, we measure time t in Eq. 
([2Q|) in units of uJq^ and introduce the dimensionless chirp 
rate a = a/uj^. With this rescaling, the initial Wigner 
distribution becomes /o = (27r)"^ exp[-(x^+II^)/2]. 
We solved Eq. ([2Q|) numerically with the same parame- 
ters as in the Schrodinger simulations and show the re- 
sults in Figs. [l][3] (subplots e-h) at the same times for 
comparison. For a better representation of the Wigner 
distributions for different nonlinear it ies, we rescaled the 
u^x axis in the Figures to v = \^u^ and ^ = \/^x. 
The dashed lines in the Figures are the separatrices, en- 
closing all bounded classical trajectories in phase space. 
We started all these simulations in the ground state, i.e. 
7 = 2, at the initial time tq = —8. Figure [T] compares 
the dynamics in phase space to that in the energy ba- 
sis in the quantum LC regime (P2 = 8), using the pa- 
rameters a 6.25 X 10-'^, ^ = 0.0042, and s 0.013. 
The pattern seen near the origin in Fig. [1] is due to the 
quantum interference with a finite number of states in 
the nonresonant region. Figure [2] shows the intermedi- 
ate (P2 = 1) case for parameters a = 10~^, (3 = 0.0067, 
and e = 0.02. Finally, Fig. [3] corresponds to the classi- 
cal AR case (P2 = 0.2) and the parameters a = 10~^, 
P = 0.0013, and e = 0.038. As weh known [21], in the 
near classical case the Wigner function becomes oscilla- 
tory on increasingly fast phase space scales. However, if 
coarse-grained (due to a finite numerical accuracy in our 
case), the Wigner function becomes almost everywhere 
positive as one approaches the classical distribution func- 
tion, despite the initial quantum-mechanical ground state 
used in the simulations. The evolution of the Wigner 
function in the last example is nearly classical with the 
quantum signature entering only via the effective tem- 
perature of the initial state. In the classical for- 
mula (p!4|) for the transition width, T appears due to 
integration over the classical Maxwell-Boltzman distri- 
bution function (see [l5j). Therefore, for the quantum- 
mechanical initial conditions, we should integrate over 
the Wigner function in a thermal state instead over the 
classical distribution. But these two distributions have 
the same functional shape, except that T is replaced by 
Teff in Eq. ([T9|) . Therefore, as also confirmed in ex- 
periments [11], one can use the classical formula for the 
threshold of the phase-locking transition at low temper- 
atures, when starting from quantum- mechanical initial 
conditions. 



B. The dynamics in the rotating frame 



di di_ dVdl ^ f^d^ 

dt^'^dx dxdu 4 du^' ^ ^ 

where, x = x/L^ u = u/ujqL^ L = \/^BT^^/muj^ , 7 = 

V = ^x'^ - ^'^x^ ^ ex COS ifd: (21) 



Here we further expand our discussion of the classi- 
cal AR limit in our system via the Wigner representa- 
tion in the rotating frame. The transformation to the 
rotating frame is accomplished using unitary transfor- 
mation (see [22]) U = ex.-p{—id^d(pd)^ where the operator 
d = {2mJiwo)~ ^ {mujox -\- ip) and cpd = J oj^dt is the driv- 
ing phase [see Eq. ([1])]. Then, by neglecting rapidly 
oscillating terms, the Hamiltonian ([1]) is transformed to 
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where 



A 



G = ^{Q^ + P')-\{Q' + P^f + ,,Q. (23) 



(22) 



The parameter /i 



3(3 



hPi\/P2 in the 



32u;o mQ;3/4 2' 

last equation is familiar from the theory of the classi- 



cal AR MJ, while A 



3/1/3 



7P2 is the dimensionless 



Plank constant, entering the commutation relation for 
the rescaled variables 



(24) 



Here Q L-^{x cos cpd + sin ^pd), P 

L~^{:^^ cos (pd — xs'mcpd)^ where = h/{muJoX), and 
the dimensionless time associated with the dynamics gov- 
erned by Hamiltonian (|23|) is r = ^/at. 

Next, we write the quantum Liouville equation in the 
rotating frame (see Ref. [S^ for similar developments for 
a constant frequency drive) 



a/ 



dGdf_ 
dPdQ 



dQdP 4 •'' 



(25) 



where D = (q^ - P^) + The initial 

Wigner distribution (p!9|) in the new variables is 



/o(Q,P) 



1 



27ra2 



(26) 



where ct^ = = ^ coth( „,^"° ). The left hand side 

of the Eq. (j25j) is identical to the Vlasov equation describ- 
ing the evolution of a classical distribution of particles 
governed by Hamiltonian ([23]) without collisions and self- 
fields. Hence, as in the fixed frame, after coarse-graining 
the fast phase space oscillations of / in the limit A ^ 
(P2 ^ 1)5 the dynamics in phase space can be treated 
classically [2l|. Therefore, both the threshold and the 
width of the autoresonant transition can be calculated 
from the classical theory as illustrated in Figs. H] and [5l 
respectively, despite the quantum-mechanical initial con- 
ditions in the problem. In other words, P2 is the measure 
of the classicality of the phase-locking transition in our 
chirped oscillator. Furthermore, in the limit of P2 <C 1, 
only two parameters, (i = |PiP2^^^ and Teff (via the 
initial conditions) fully characterize the AR transition. 
This result is in agreement with Eqs. (p!3|) and (p!5]) for 
the AR threshold and its width, where, remarkably, ja 
and Teff enter separately. 



VI. CONCLUSIONS 



In conclusion. 



(a) We have studied the interrelation between the 
quantum-mechanical and classical dynamics of phase- 
locking transition in a Duffing oscillator driven by a 
chirped frequency oscillation. We studied the conditions 
for a continuous phase-locking in the driven system, such 
that the energy of the oscillator grows to stay in reso- 
nance with the varying driving frequency. The problem 
was defined by the temperature T and three parameters, 
i.e. the driving amplitude e, the driving frequency chirp 
rate a, and the parameter P characterizing the nonlinear- 
ity of the oscillator. The nonlinearity in the problem was 
essential, since no persistent phase-locking in the system 
could be achieved for 13 = 0. 

(b) We have exploited a more natural representation of 
both the quantum-mechanical and classical dynamics in 
the problem via just two dimensionless parameters 

Pi = e/\^2mhuJoa and P2 = {3h/3) / {4:m^/a) , instead of 
£, a, and p. We have shown that P2 describes the classi- 
cality of the phase-locking transition in the system, such 
that, for P2 <C 1, the system arrives at its classical au- 
toresonant (AR) state after passage through linear res- 
onance even when starting in the quantum-mechanical 
ground state. In contrast, for P2 Pi + 1, the transi- 
tion involves the energy ladder climbing (LC) process, i.e. 
a continuing sequence of separated Landau-Zener transi- 
tions between neighboring energy levels. The parameters 
Pi ^2 have a meaning only in the case of a finite chirp rate, 
which introduces a new time scale, Ts = l/v^^ the 
problem. 

(c) The probability of transition to the phase-locked 
state versus Pi has a characteristic 5'-shape (a smoothed 
step function). The value of Pi yielding 50% transition 
probability can be viewed as the threshold for the phase- 
locking transition. We have calculated this threshold and 
its width in both the quantum-mechanical LC and clas- 
sical AR limits and compared the results to those from 
quantum-mechanical calculations starting in the ground 
state of the oscillator (see Figs. |4]and[5j). We have found 
that, while in the LC limit the threshold is independent of 

P2, in the classical AR regime, the threshold is defined by 
1 1/2 

the combination ji = ^PiP2 of parameters. The agree- 
ment of the theory and simulations in both limits was 
excellent, but characteristic oscillations of the threshold 
and the width were observed in the intermediate regime 
1< P2 < 5. 

(d) We have also studied the dynamics of the phase- 
locking transition in phase space by using the Wigner 
function representation, to explain the quantum satura- 
tion of the width of the threshold for AR transitions. The 
analysis of the Wigner (quantum Liouville) equation in 
the chirped rotating frame clarifies the role of P2 as char- 
acterizing the degree of classicality in the phase-locking 
transition problem. 

(d) A possibility of engineering and control of a desired 
quantum state of the oscillator via ladder climbing pro- 
cess (see an example in Fig. [1]) seems to be attractive in 
such applications as quantum computing. A generaliza- 
tion of this study to include possible decay, decoherence. 
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and tunneling processes also seems to be important in 
future studies. 
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